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Abstract- The present paper develops a model for calculating the 
rate of photonuclear reaction (beta decay) for nuclear 
transmutation in the reducing of radioactive waste. A 
photonuclear reaction is viewed as an incident photon creating 
superconducting hot spot (hot belt) across a nucleon from the 
composition of an unstable nuclide (radioisotope), followed by a 
thermally induced vortex crossing which turns superconducting 
hot belt into the normal state (vacuum) resulting in a vortex 
assisted photon beta decay. Because of this model requires data 
on energies of vortex and on currents from inside of the nucleon, 
an analog of a superconductor model was developed for the 
nucleon. Thus, it revisited the dual Ginzburg-Landau model for 
the calculation of Lorenz force, monopoles current, and the 
energy of vortex lines for a vortex triangular lattice type 
Abrikosov within a nucleon to find their meaning. For now, it 
was found that these energies would correspond to the subatomic 

particles, w,z, Higgs bosons, pion 71* , and nucleon itself. As 
check points of this model are considered, the deduction of the 
fusion temperature of two nucleons, even this model itself for 
photonuclear reaction mechanism, the explanation of the 
mechanism of natural beta decay as dark counts and comparison 
with the results of the IBA nuclide substructure model are also 
considered. Finally, the model provides counts (decay) rates 

( R ), and allows the estimation of intensity « ( photons s), and of 
threshold value of photon energy (MeV), in order to 
achieve R Ir =1 , or ( 100% ) efficiency. Such efficient 

pc I h 

installations could be, for example, the laser ELI-NP in 
construction in Romania for the high energy photon source. 

Keywords-Beta Decay; Nuclear Transmutation; Photonuclear 
Reactions; High Energy Lasers; W,Z,H Bosons; G-L Theory 

I. INTRODUCTION 

The Photonuclear reactions are proven to be of potential 
interest for nuclear transmutation [1] in the handling of 
radioactive waste, which is a big problem in nuclear energy. 

Photonuclear cross sections are quite large (~ 1 barn) at the 
giant resonance (GR) region of E ? =15H-30MeV • 
Major (y,2n) and (y,n) reactions on long-lived nuclei may 
transmute them to short-lived or stable nuclei. 

Till now, the coherent photonuclear isotope transmutation 
(CPIT) produces exclusively radioactive isotopes (RIs) by 
coherent photonuclear (y,n) and (2y,n) reactions via El 
giant resonances [1]. Photons to be used are medium 
energy E(j)=l2-25[MeV] photons produced by laser photons 
backscattered off GeV electrons. 

The photon intensity of around 10 8 12 Is is realistic by using 
intense lasers and intense electrons. In future high intensity 

photon sources with the intensity of the order of io"" 5 /i will 
be possible. Therefore, there are two conditions for a 



photonuclear reaction to happen: the intensity, and the 
threshold energy of photons, respectively. 

Otherwise, in [2a] in the frame of nanowire single -photon 
detector (SNSPD) research, the transition from a current- 
biased metastable superconducting ( s ) state to the normal ( n) 
state is triggered by vortices crossing the thin film 
superconducting strip from one edge to another due to the 
Lorentz force. Thus, a photon creates a spot which suppresses 
the order parameter that lowering the energy barrier for vortex 
crossing. 

The model developed in [2a] is based entirely on the 
Ginsburg-Landau theory of superconductivity. 

In the present paper, based on the above models [1, 2a], a 
photonuclear reaction is viewed as an incident photon creating 
superconducting hot spot (hot belt) across nucleon from the 
composition of an unstable nuclide (radioisotope), followed 
by a thermally induced vortex crossing, which turns 
superconducting hot belt into the normal state (vacuum) 
resulting in a vortexassisted photon beta decay. 

Since a model requires data on the energies of the vortex 
and on the currents from inside of the nucleon, an analog of a 
superconductor model was developed for the nucleon. 

Usually, the masses of W,Z are calculated by taking into 
account and a priori Higgs field, and the default using the 
Higgs mechanism and Higgs boson. 

Soon after the advent of QCD, Hooft and Mandelstam [3] 
proposed the dual superconductor scenario of confinement; 
the QCD vacuum is thought to behave analogously to an 
electrodynamic superconductor but with the roles of electric 
and magnetic fields being interchanged: a condensate of 
magnetic monopoles expels electric fields from the vacuum If 
one now puts electric charge and anti-charge into this medium, 
the electric flux that forms between them will be squeezed 
into a thin, eventually string-like, Abrikosov-Nielsen-Oleson 
(ANO) vortex which results in linear confinement. 

The dual superconductor mechanism [3] is an alternative 
that does not require the ad hoc introduction of a Higgs field 
but instead uses dynamically generated topological excitations 
to provide the screening supercurrents . For example, U (1) 
lattice gauge theory contains Dirac magnetic monopoles in 
addition to photons. The dual superconductor hypothesis 
postulates that these monopoles provide the circulating color 
magnetic currents that constrain the color electric flux lines 
into narrow flux tubes. Hooft has shown that objects similar to 
the Dirac monopoles in U (1) gauge theory can also be found 
in non-Abelian SU (N) models. 

The results are consistent with a dual version of the 
Ginzburg-Landau model of superconductivity. What is 
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important in understanding field (magnetic) dependence was 
Abrikosov's field theoretical approaches based on Ginsburg- 
Landau theory [4] for type I superconductors 

( K > k—XlS, , A is the penetration depth, £ the 

coherence length ) and type ones ( k >l/V2 ) II 

superconductors, which allows magnetic flux ® to penetrate 
the superconductor in a regular array quantized in units of 
elementary flux quantum ($>=7rticle .The importance was the 

n+ , (n=0,±l, — ) 

V 2) e 



quantization in a ring , flux 0= 



In the present paper we revisited G-L model [4, 5, 6, 7, 8, 
9, 10, 11, 12, 13] in order to calculate the values of the 
Lorenz force, the current, and the energies of the Abrikosov 
vortex lines inside of the nucleon, in natural units, firstly, in 
view to search for its relevance; for the time being, it was 
found that this would correspond to energies for subatomic 
particles, such as that of W,Z , Higgs bosons, and of pion 

, and secondarily to be used in the model of vortex- 
assisted single photon count. 

In Section IV the model is applied to formulate the 
mechanismof natural beta decay as dark count. 

In this model to a superconductor analogue, we do not use 
any a- priori field (Higgs), we use only a comparison with the 
electrical field generated by the pair (dipole) qq . 

If these kinds of predictions continue to be verified by 
appropriate experiments, one may have a cheap solution to 
remove the radioactive waste of used-up nuclear fuel rods of 
fission reactors in a reasonable time period. 

II. DESCRIPTION OFTHE ANALOGUE MODEL OF NUCLEONT O A 
SUPERCONDUCTOR 
The normal cores that exist in type-JJ superconductors in 
the mixed state are not sharply delineated. The value of 
number density of superelectrons n is zero at the center of 
the cores and rises over a characteristic distance £ , the 
coherence length. The magnetic field associated with each 
normal core is spread over a region with a diameter of 2A, 
and each normal core is surrounded by a vortex of circulating 
current. 

The QCD vacuum can be viewed as a dual superconductor 
characterized by a monopole condensate [3, 9, 11], when 
embedding a static qq pair into the vacuum. The core of the 
flux tube is just a normal conducting vortex which is 
stabilized by solenoidal magnetic supercurrents, j , in the 
surrounding vacuum. 

In order to calculate distinctly the energy states (masses) 
in natural units, firstly we re -derive the field equations of 
magnetic monopoles current, of the electric flux and energy 
states. 

Therefore, here is adopted a basic dual (B <=> E) form of 
Gin zburg- Landau (G-L) theory [4, 5, 6, 7, 8], which 
generalizes the London theory to allow the magnitude of the 
condensate density to vary in space. As before, the 
superconducting order parameter is a complex function ^(j) , 

where [^(i)] 2 is the condensate density n . Also the wave 



function is defined as y/(x) = yfn^ exp(i0(x)) , where „ is the 
London (bulk) condensate density, and ^are real functions 
describing the spatial variation of the condensate. 

The characteristic scale over which the condensate 
density varies is the G-L coherence length or the vortex 

core dimension. The x denote the radial distance of points 
from the z -axis, the superconductor occupying the half space 
x > ■ Outside of the superconductor in the half space x < o, 
one has b = e = h = H > where, "the external" vector H is 

parallel to the surface. The >ytheory of superconductivity [3], 
[4] is an application of the Landau theory of phase transitions 
to superconductivity. In this case, some scalar complex v 
function fulfils the role of the order parameter. 

The final expression for the free energy then takes the 
form. 



/=/»+! 



h 1 



2m 



-i — A \y/ 
he 



, ,2 b, ,4 B 
a\w\ + — \w\ H 

11 2 1 1 Sn 



>dV 



(1) 



One can obtain the basic equations of Gin zburg -Landau 

theory by varying this functional with respect to A and y/* . 

Carrying first variation with respect to A, we find a simple 
calculation: 



Sf=\ 



ieh 

c (if/ V ' tp-yN y/ ) - 

2m 

2e 2 , ,2 curlB 

\y/\ Ah 

m An 



SAdV+{div(SAxB)— = 
J V ' Att 



(2) 



To minimize the free energy, the expression in the 
brackets must be equal to zero. This results in the Maxwell 
equation. 



An 1 
curlB = j s = — — j s (in _SI) 

c c s' 



or 



An 

Vx Vx A = — .j s 

c ' C'E, 



1 

2— J. 



provided that the current density is given by 



ieh, . ,\ 2e , , 2 . 

— I y/ V y/ - yN y/ H \y/\ .A 

2m m.c 



(3) 



(4) 



(5) 



The following Equations (6) and (7) form the complete 
system of the Ginzburg-Landau (G-L) theory. 



An 1 

VxVxA = ./ =——j 

J x 2 J s 

C C £ „ 



(6) 



According to the definition of we use y/ = yj~n~ exp(i^) . 

The second integral is over the surface of the sample. 
The volume integral vanishes when 

h 2 ( 2e Y , ,2 n > 

V-i — A y/ + ay/ + b\y/\ y/ = (') 

Am { he J 
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Where ^ _ s o- m - c = i.i7e — 16[m] > an( ^ tne quantized flux: 

v n ,.e J 

<j) = ^£,and |^| 2 =n s ; « = 3 _monopoles/v *l.e-45m 3 , 

e 

V =4/3;zr 3 = 2 r = 0.48[./m] f = 8.8e-12[C 2 JV-'.m 2 ] 

Since the magnetic charge [17] of monopole being 
tic AneJic 137 



U =4^o 



2e 



68.5e 



and assuming that the classical electron radius be equal to "the 
classical monopole radius" from which one has the monopole 

2 = 4700m , and the value of X remains 



mass m M = g]mj 
unmodified. 



The superconductors of second kind are those 
with K > l/V2 , and X > £ . 

We now consider the phase transition in superconductors 
of the second kind. 

For this we can omit the non-linear (|^| 2 y/) term in (7), 
1 



we have 



Am 



2e 

-iti A 



¥ = \ a W 



(8) 



This equation coincides with the Schrodinger equation for 
a particle of mass 2m and charge 2e (in the case of dual, the 
factor 2 for the charge, which is specific to the "pairs", it is 
actually 1) in a magnetic field H Q (m our case the chromo- 

electrical flux £"(0) )■ The quantity |a| plays the role of 
energy ( Ey/) of that equation. The minimum energy for such 
a particle in a uniform electro-magnetic field 



H n is an "external" 



1 1 2eH n 
is e m =-ha> a =-tt- =>[/] 

2 2 8mc 

electro -magnetic field of a dipole created by the pair uu (the 
chromo electrical colours field). 



H=E = 



de 



Ansj 



■ = 8.33e24 



TV 

c 



(9) 



where r = 0.05[/m] is the electrical flux tube radius, 
d — 0.48[/m] the distance between the two quarks 
charges, usually H[A/m] , but here is used as 

B = uH 



Am 



Hence, Equation (8) has 
|a| > 2h *eH /8mc 



a solution only if 



or 



H <- 



TltlC 



Amc\a\ 



he 



■<H_ 



4> 



2< 2 



(10) 



- = 8.33e24 



and in terms of B , = H „ = 



nfic 



= 2.7el6 



particle energy is X > £ with £ 



2^ 2 .c 

X 0.117 



.The 



K 



1.05 



Am 
0.1114[ym], 



or k > 1/V2" > 1 = 1.05 (of type H-superconductor). 

One of the characteristic lengths for the description of 
superconductors is called the coherence length. For example, 
a transition from the superconducting state to a normal state 
will have a transition layer of finite thickness which is related 
to the coherence length. 

More exactly, this quantity is called the correlation or 
healing length [4], and is defined as 



1/2 



\T-Tj 



(11) 



where £ = atio F j ' E g is for T — > , a — 2/k from [19], 
E -gap energy, k Boltzmann constant; at confinement 



T = 175[MeV] 2el2[K] 
electrons (monopoles) is 



o r =■ 



and the Fermi velocity of 
^2*A10QmE F 



4700m 



(12) 



where as the Fermi energy we have for monopoles condensate 
viewed as bosons condensate 



h 2 



n n 1 



4700m 



= Tk B = 0.1E F 



where 



E f = 

2* 4700m 

( 

numerically, 



{bnn s ) 



2/3 



(13) 



1.06e-34 
2 * 4700 * 9.e - 3 1 



(3ff 2 *3/(V*l.e-45)) 



we have: 
: 9.32£--12[/]^55[M«V] 



where V = 1.5[/m]\ and the velocity of monopole is 

u F = 0.55e8 < c = 2.997 eS[m/s] 

and 



^ =0.6*0.7- 



.1.06c- 34 *0.55e8 
1.38e-23*2.el2 



:1.02e-16[m] at r (14) 



Usually, the flux quantum is defined as: 

O =nhc/e-> usually — = 2.01e-l5[Tm~ 

e 



(15) 



In [4], the G-L equations are solved analytically only for 
X » £, (near T this means k » 1 ). 



Thus, in the range ' 



(16) 
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The induction B is given by: B(x) 



2<S 



^log 



V X 



(17) 



This equation is valid only at all distances x » £, (18) 
The solution at x — > oo is B(r) = const -K (x/ A) , 
where K is the Hankel function of imaginary argument. The 
coefficient must be defined by matching with the solution of 
(17). Using the asymptotic formula K (x) « \og{2jyx) 

for x « 1 , where y — e = 1.78 ( C is Euler's constant), we 
20„ 



finally have B(x) = 



-K (x/A) 



2nA c 

Using equation (19), we can rewrite (17) as: 

2A 



B(x) 



20„ 



-log- 



x « A 



(19) 



(19.1) 



2kA c yx 

In opposite limit of large distances one can use the 
asymptotic expression K (x) « (nfl x) 1 ^ 2 e~' x forjc>>l. 

Thus, at large distances from the axis of the vortex line the 
field decreases according to 



B(x) = 



20„ 



c($7TXA ) 



3 s 1/2 



, x » A 



(20) 



Accordingly the superconductive current density decreases 
(in SI): 



c dB , „ x 2c 7 s n <I> n 

-{47TC£ ) = 



4-7T dx 



8(2^- 3 xA 5 y 2 c 



-x/A 



(21) 



We can now calculate the energy s of the vortex line; we 
obtain for the energy per unit length of vortex line. 



■■ c e n 



( 2 <0 








log 




\\nlc ) 







Here, the factor 
from (cgs) -> (SI) . 



\ncs n 



(22) 



used to convert 



Because the magnetic induction of the monopoles current 
which is powered by electric field given by a pair of quarks 

( H ), it has the raw flow consequences squeezing this 

cromoelectrical flux into a vortex line, followed by forcing an 
organization into a triangular Abrikosov lattice, see Figs, la, b. 




Fig. la Abrikosov 's triangular lattice for anucleon (proposal) 




Vortex (W) 



Giant- Vortex 



Fig. lb The Giant-Vortex (after ref. [24]) type arrangement for the nucleon 
(only as an illustration) 

The core of every vortex can be considered to contain a 
vortex line, and every particle in the vortex can be considered 
to be circulating around the vortex line. 

The presence of vortex line increases Z, Vortex(W), the free 
energy of the superconducting media with sL , thus it is 
thermodynamically favorable if the contribution is negative; 

i.e. if sL-<3> /c-H L/4xjU<iO , and B =— 1 , 

Mo 



/Z =l/c 2 £ Q 



or H n < H . —■ 

cl 



(23) 



Substituting (22) in (23), we find the lower critical field 



B=H 



2^>l 2 



-log 



nfic 

— r log(*r) = l.el5 

TtX C 



J 



Am 
(24) 



E, — 0.1 1 14 , and when near the axis, for x = 0.116 = <^ 

from (20) when the induction is B(%) = 2.e\5 = 2H cl (25) 

The energy of interaction of two vortex lines is separated 
by a distance d from each other is given as: 



20 



40 



B,(x)- 



K, 



(26) 



\n " %nA l 

One can use the asymptotic expression for £. (see Eq. 
(27)). 

401 



s . = 



~7/2 3/2,2 

2 ' TV 1 A 



d J 



, x»A (27) 



When the distanced w A, » E, , the cores of vortex lines 
overlap [4]. 

Let us consider a closed contour near the surface of the 
cylinder. The change of wave function on passing round the 
contour is 2kvS , where S is the cross-section area of the 
cylinderand v is the number of vortex lines. The electric flux 
is 

2m r j, 







20 o v5 



eh 



dl 



(28) 



Let us recall that similarly relationship [4, 5]. It was 
introduced for the first time by London, called fluxoid 
equation. 
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Each fluixoid, or vortex, is associated with a single 
quantum of flux represented as <D , and is surrounded by a 

circulating s uppercurrent , j , of spatial extent, A . As the 

applied field increases, the flu xo ids begin to interact and as 
the consequence ensembles themselves into a lattice. A simple 
geometrical argument for the spacing, d of a triangular 
lattice then gives the flux quantization condition [13], 

S/=^<E (29) 
V3 

where B is the induction. 

The solution of Ginzburg-Landau phenomenological free 
energy [13] is useful for understanding the Abrikosov flux 
lattice. The coordinate-dependent order parameter describes 
the flux vortices of periodicity of a triangular lattice. 
Fluctuations from tj) change the state to^/- ; the minimization 
of free energy with respect to iff , gives the ground 

state <f>(r/0) . 

The free energy is given by, 



f = l + \ 



h 1 




2e 




(• 


-i — A 


2m 




tic 






b, |4 


+a | 







(B-H,y 

87T 



>dV 



(30) 



The average magnetic induction is B(— y,0, 0) ■ The free 
energy has solutions of vortices of triangular form The 
coordinates of the three vertices of a triangular vortex are 



given by(0,0),(/,0), and 



2 2 



/. The fluctuation from 



ground state corresponding to that of triangular lattice is that 
for small fluctuations. The deviation of the free energy from 

the mean-field value F — F FM with respect to the thermal 

energy, k B T , can be used to obtain the physical properties of 

the fluctuations which are useful for understanding the melted 
vortex lattices. The deviation from the triangular Abrikosov 
lattice is defined as 



D = <|^-fl 1 ^(r/0)| 2 )/fl 1 2 



(31) 



which uses the spatial and thermal averages calculated with 
the probability exp (— F/k g T} . Classically, 



D 



k B T 



F-F„ 



(32) 



measures the fluctuations fromthe triangular vortex state. The 
fluctuation in the distance between vortexes becomes: 



-casei, (i-r FM /r) = io-V 4 V/ 3 



(33) 



-Case 2, l-T FM /t = c~ l B 

7 FM I c 



5/4 



(34) 



-Case 3, a vortex transition below the transition temperature 
see [13], 

where, T FM , the flux-lattice melting temperature, and 
c = 0.1 from Lindemann criterion of lattice melting when 
d 2 —c 2 l 2 , the flux quantization condition Z 2 =<D /B , 

, ey[2 2 

B = 2n;n/>c , where k = A'H(T) , and 

he 

H = Kiln k.H , . 

For numerical values T = \15[MeV] , in case of 
symmetry breaking, in case 1, it results T FM = T c , and in Case 
2, it results T FM = 100[feV] by using (29) in place of 
O /fi = d 2 with d = 0.3982[/m] (a very precisely value), 

and tc = 1 , which is the temperature of fusion (melting!) of 
two protons . 

This triangular lattice corresponds to the arrangement of 

the quarks pairs uu,uu,dd in the frame of a nucleon, see 
Fig. la, b. 

A direct numerical analysis allows obtaining the following 
values for the current, force and energy. Thus, from (21) the 
current is given by: 



7 ;=1.15e7[A/> 2 ] 



(35) 



where x = X . 



For x = 1- X , the current density decreases at 
j F =3.e5[A/fin 2 ]- 

Note that velocity v F , moreover, if one considers the 

monopole current given by equation as j = n v F g D , where 
the magnetic charge is : 

tic _ AkeJxc 137 
2e~~ 



- e = ^ e = 6&.5e (36) 
2 



If we use the range x = 0.l<A , then the current is 
obtained by derivation of (19. 1): 



7<* — 



1 1 1 \ntic c 2 £{) _ 1 1 1 
4xA 2 e c ~4xA 2 



V 



^c = ^v Fi g D ^> (37) 



^IL = — = 0.302e8[m/s] 
c 9.89 ' 

The Lorentz' force to squeeze the flux tube is: 

F L = qv Fi B = 1.6e- 19 * 0.1 * 2.992e8 * 4.7el5 = 2.25.e4[Af] 

(38) 

when B is given by (19.1) and x = A , for the upper limit: 



B(x) = 235el5[j/Am 2 ] 
With B from (20) and for x » A , we have 



(39) 
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B{x) = 3.8e-16 



Am' 



(40) 



where x — 12X — 8.4[/m] , then, the force becomes 
F L = 1.8e - 26[N] , or in terms of energy 



e = F * x = 945[MeV] 

barrier L L J 

or the nucleon overall. In case of x = E, — > (0) , in ( 17), 

./ 



B(0)=1.03el5 



Am 



which respect (25). 



The magnetic energy results from (22), and 
(A » x »£) from (16): 



(41) 



(42) 



for 



£ = 1.09e-10[7/>i]=>0.66[GeV/>i] (43) 

The force on the flux tube which corresponds with the 
well known "string tension" defined in refs. [9-1 1]. 

Now, from (27) and d « (4 h- 6)/l »^,we have 

f im = 2.3e - 1 l[j/fin] => lU[MeV] (44) 

What would be the value of the mass of the pion n + , 
composed of a pair of quarks ud interacting at a distance 
d ~ 2/3 = 5.65* X = 0.66[fm] of the radius of the nucleus. 

Now, others important values of energy, as will be 
necessarily in following, it could be calculated: 



s (0) = s iBt (d = x-A; 
x = 0. 14) * 0. M[fin] = \.e- 09[/] 
and from(38) with x = <!; =0.101 [fin] 

s Qh =Vc 2 s () (2H cl ) 2 /8x = 5.e-n[J] 
Now, the vortex energy is : 

^=VcVC/8*- = 1.16e-08[7] 



(45) 
(46) 

(47) 
(47.1) 



where V is the volume, see Fig. la, b. Accordingly, the 
corresponding equivalently masses are 

M = s mmx jc => 73[GeV] , which seems to be equal to the 

mass of W boson. 

The energy of the neutral boson Z is assimilated with the 
vortex-vortex three pairs interaction energy [24], 

£ =3*s. =9l[GeV], when from (27) 

Z mt- pair L J ' v * 



e^-^eJA = 0.117-0.116; 



(48) 



x = 1.42) = 4.85e - 09[7] -» 30.33[GeV] 

is the energy of each of three pairs of vortex outermost 
(d = 0) vortices lines which interacting (repel) at the center 

of the triangle situated at x = 1.4/1, thus, being generated a 
neutral current in the zone of Z during the triangular 
arrangement of the lattice, see Fig. la, b. 



Now, it is possible that the vortices start to coalesce into a 
giant vortex (GV) [26], see, Fig. lb. , 

Thus, from (27), results another energy state-maximum 
possible ( d = ), probable that of Higgs boson (H): 

s H = 3 * e u (d = 0.117 - 0.1 16; x = X) = 2.17 'e - 08[7] -> 135[GeV] 
Notice that this is not in fact a particle, since it contain 
others subparticles (W ,Z,u,d,monopoles ), so during high 
energy protons collision (CERN), it cannot be obtained as 
itself. 

Here, a factor of 2 was introduced to correct on 2e for 
"pairs" in the G-L model. 

JJI. THEPHOTONUCLEAR REACTION MECHANISM IN G-L THEORY 
As in [2a], we derive the energy barrier for vortexcrossing 
and crossing rate from the standard Ginzburg-Landau (GL) 
functional with respect to the superconducting order 
parameter y/(r) (normalized to its zero-field value in the 

absence of current) and the vector potential A in the presence 
of monopoles current induced magnetic field H ■ The Gibbs 
energy is given by (30). 

Like in [2], we assume that in sufficiently narrow 
strips, w = r , (a) the cloud covers the entire width of strip 

when it reaches its maximum size; (b) the quasiparticle 
density in the cloud (hot belt) is close to uniform; and (c) 
quasiparticles suppress the superconducting order parameter 
inside the hot belt, but their density is not sufficient to convert 
the hot belt to the normal state. Thus the superconducting 

condensation energy density H , in the hot belt satisfies 
the inequalities < H^I&TT < H 2 c(Q J%7r . 

It follows that the vortex crossing rate via hot belt 
(denoted by subscript h), R uh {I), is enhanced in comparison 
with that for dark counts, because the parameter 
v = s Ik T » H 2 link T is reduced in comparison with 

h Oh I B ch 1 B L 



v =£ /7\and s o 



cD / 



e = l.e-09[7],from(45). 



There are three superconducting current-biased regions to 
be considered, which are separated by three characteristic 

currents I* < I , < I „ . Currents below /* counts are absent, 

h ch cO h 

because vortex crossings do not result in the formation of 
normal-state belt. The "hot" current /, is determined through 

h 

the energy balance for destroying the superconducting 
condensate in the hot belt similar to that for dark counts, 

however, with the thermodynamic field H < H (0) = H i , 

from (24). 

For lower currents in the interval / >/>/*, the rate of 

Ch ii 

photon counts shows a power-law current dependence, 

R (/)«/"'*' 

Here, we model the temperature dependent parameters, 
I ch (T) = I ch (0)s 3 T , v i (T) = v i (0)s 2 T /T , with 



e T =(l-T/T e ) 



1/2 
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For currents / > I c /3, the energy released by phase slip is 
sufficient to destroy superconductivity. 

In the region /*</</* = I c0 /3, vortexcrossings via the 

hot belt destroy the superconducting state inside the belt 
during the time r of the hot spot's lifetime. During this time, 
the belt is in the normal state, which causes the "vortex- 
assisted" photon count. 

A sketch of the strip and of the belt across is shown in Fig. 

2. 



Pie ) 




Hot spot 






Fig. 2 The proposed photonuclear mechanism: from left to right, illustration of 
incident photon creating superconducting hot spot (hot belt) across nucleon, 
followed by athermally induced vortex crossing together with an electron (bias 
current), which turns superconducting hot belt intothe normal state resulting in a 
vortex assisted photon beta decay 

If />/>/, the rate of single photon counts, R , is 
the same as the rate for hot spot formation r , because the 
barrier for vortex crossing is now absent. For hot counts R , 

we have R = rjR , where R is the photon rate and rj t is the 

quantum efficiency of hot spot formation caused by a single 
photon. 

The vortex assisted photon count rate from [2a] is: 
R pc =R h [l-exp(-K)] 

where: 



(49) 



Ak B TcR £ 



7Z®„ 



V 






v h +1 


\2n ) 









9? 



=i co{ v h/ v o f 2 



(50) 

(51) 
(52) 



In the case of ft disintegration n — » p + e +V , 
or udd -» uud + e~ + v e + W (S0GeV) , 
or fi?(-l/3e) + 2*3/3e = M(+2/3e) + e, 

or e + u(2/3e)-2*3/3e = d(-\/3e) (53) 

When the vortex equilibrium is disturbed by the 
transformation of one quark ( d ^ u ), it will be accompanied 
by a release of a W boson or the crossing vortex in this model. 

We can suppose than along the hot belt induced by the 
incident photon, the charge 2e + W creates a bias current 



( 7>2/3e[(l/(v, i /v )' /2 ) = 3]=>2e ) who circulates due of the 
potential difference between the vortex and the rest of isotope. 

At the first sight, the ohmic resistance of this ad -hoc 
electrical circuit created by the bias current is given as: 

R = -?- (54) 

r V 2 

GL vertex 

where the vortex potential is V =H^, Z/ o from (30.1), 
and the power is U „ T = S t IT , where U „ is given 

r PI vonexiW 1 )/ P ° 

by (47.1) , and t cl = xfr/(MT) = 1.5e-24[s], the Ginzburg- 
Landau life time of M /± bosons. 

Numerically, with T = \15MeV , E =k T , result in 
v =s w /kT =l.e-09/(1.38e-23*2.el2) = 36.2 

^from (43); Vh =S ^>^l E ph =5-g ~ 1 1 / E r'»° • 
from (46). 



where 
i' rag , where ^ int 



The value of E is determined by trials in order to 
have/? IR. =1. 

pc I h 

Since, a photonuclear reaction could be viewed as a 
vortex assisted photon beta decay, and with the above data 
related on energies of vortex and on currents from inside of 
the nucleon, are calculated the beta decay rates. 

Now, let us discuss the m Mo — > mm Tc isotopes, which are 
widely used as SPECT isotopes [1]. The 

" Mo (ft~,T,= 65 h,Q = 1.356MeV ) isotopes are produced by 
m Mo(y,n) reactions and also m Mo(y,p) reactions followed 
by /T decays. 

From our model results, only from B~ decay chain 

( *Mo -> M ™7c ), a decay rate of R pc = 4.7 x lfj 7 counts/ s 
is attained, with the incident photon rate of 
R : = I. el3 photons /s and of ~ l5MeV photons energy 

(or r Ir s 4.7 x io '). Therefore, in order to transmute entirely 

( R I R : = 1) it needs, or an irradiation of 60h = Ti = 65h , 
either, to proceed much faster (instantly), the photons energy 
must be of a very precisely threshold value > E = 20Mev , 
see Fig. 3. 



Rpc/Rh-*-RV-Rpc 




Fig. 3 The vortex -assisted photon count rate R j R vs. bias current given 
inEqs. (49 - 53) 
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137 .r, v 137 D , - , — 

5S Cs — > 56 Ba + e +v 



These are the main results of this model, namely, a 
threshold value of photon energy, for an incident 

rate R = R = 1 x 10" [photons / s] , or in other words, the 
necessary parameters of photons source. 

Notice that, at these parameters, the same type of instant 
rates it happens for all beta-decay isotopes, like: 85 A> , 
Q p = 0.687MeV and T ]/2 = 4Sh 

y-ray = 5 UKeV (0.46%) T l/2 = 10.756yr ; 

y - ray = 15LKeV(95%) T l/2 = 4.48/z 

Q=\.\15MeV 

P 

y-ray = 661KeV(85%) ^T i/2 = 30.8yr etc., i.e. these rates 
are not dependent of the nuclide type, but only on E . 

As it was mentioned in [1], CPIT with large efficiency is 
quite attractive from ecological view points. Since 
the GeV electrons stored in a storage ring lose little then- 
energy via interactions with laser photons, they remain in the 
ring. The laser electron photons are efficiently used for 
production of desired radioisotopes ( RIs ) by (y, n) reaction 
when a valence neutron of the nuclide is dislocated. Then the 
overall efficiency of the RI production is larger than that of 
the charged particle accelerators and nuclear reactors. Fast 
neutrons are used for nuclear transmutation but only a small 
fraction of the charged particles are used to produce the fast 
neutrons. Intense electron accelerators provide bremsstrahlung 
photons for photonuclear isotope transmutation but most of 
photons are below the threshold energy of the photonuclear 
reaction. There are several programs of intense photon 
sources under progress MAGa-ray, such as the project at 

LLNL for high intensity photons with 10 12 photons/sec in the 
MeV region and so on. ELI-NP is the one at Romania [35] 

for higher energy photons with 10" photons/sec in the GR 
energy region of E(y) > 19.5[MeV] ■ They plan to achieve the 

intensity of around \0 lS / s . In Tokai Japan, the ERL (Energy 

Recovery Liniac) project is under progress to provide intense 
photons for resonance fluorescence. 

A. Beta Decay Halftime Calculation with a G-L Model 

Below we will demonstrate that the mechanism for beta 
decay of radioisotopes is the same as the dark counts in the 
case of superconductors [2b]. 

In [2b] we derive the energy barriers for three dissipative 
processes mentioned within the GL theory. In the following, 
by analogy wit ref. [2b], we consider a thin-film strip (one of 
three vortexes of the nucleus) of width w = r * X (see Fig. la, 
b). We choose the coordinates so that it is < X < W . Since 
we are interested in bias currents which may approach 
depairing values, the suppression of the superconducting order 
parameter (y/) must be taken into account. Also as in [2b], it 
used the standard GL functional, given above in (1). 

A vortex crossing from one strip edge to the opposite one 
induces a phase slip 2ft - <p (a every phase slip meaning O 



energy released) without creating a normal region across the 
strip (one of three vortexes of nucleus) width. 

From [2], the asymptotic estimate for the dark counts rate, 
results as: 



+kJc 2 Rjm>t?(^ ( 



V 



(55) 



0/ 



and where the bias current is : 



I = —I K(l-K 2 ) 



8^A 



(56) 



(57) 



Here, h 7 = X -the axial ( z ) height of the monopole 
condensate. 

Here, the critical current at which the energy barrier 
vanishes for a single vortexcrossing: 



I = 



2ju 2 wl 



2.72^ 

And the thermodynamic critical field is: 



H. = 



2y[2x£Xc 



(58) 



(59) 



where ju 2 = 1 — K 2 , and 



Y( Z ) = (\ + z 2 f +l)l2 exp [vz tan- 1 (l/z)] (60) 

where u h = T GL (e vortex - Q bind )/tl is the energy of the 

vortexduring crossing the barrier of height s — Q. ., by 

quantum tunneling in place of the thermal activation used in 
[2b], and overpassing an ohnic resistance along a transverse 
path way of the nuclide: 



q nuclide 



R+2^/w) 2 R 



(61) 



nuclide 



Here ' £ vortex = M W from ( 41A ^ and Qbind iS the bet£l 

decay energy as obtained from the data of each radionuclide 
of beta decay type (see Chart of Nuclides 2010). 

In the case of beta disintegration n — > p + e~ + V , or 
udd — > uud + e ~ + v e +W~ (SOGeV) and the bias 
current is : d(- l/3e) + 2* 3/3e = u(+2/3e) + e. 

In /? + decay, energy is used to convert a proton into a 

neutron, while emitting a positron ( e + ) and an electron 
neutrino ( v e ): Energy + p — > n + e + + V e - 
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In all the cases where /3 + decay is allowed energetically 
(and the proton is a part of a nucleus with electron shells), it is 
accompanied by the electron capture (EC) process, when an 
atomic electron is captured by a nucleus with the emission of 

a neutrino: Energy + p+e~ —>n + V e . 

Therefore, the ad-hoc bias current created during vortex 
crossing through the energy barrier is: 

e + u(2/3e) - 2*3/3e = d(- 1/3 e) 

At the first sight, the ohmic resistance of this ad-hoc 
electrical circuit created by the bias current I(e+) due of 

quarks transformation (d —>u), or (EC(u — > d) , is given 



as: R 



nuclide T , 2 

GL vortex 



(62) 



and the superconducting quantum resistance is 
R q = n/(2e) 2 = 6.5&Q , where the vortex potential is 

V vortex =H Q ^ tf from (9). 

Giordano [21] has suggested that phase slips due to that 
macroscopic quantum tunneling may be the cause of the low 
temperature resistance tail in the ID wire he studied in zero 
fields. One possible mechanism for our low temperature 
resistivity tail could be quantum tunneling of vortices 
through the energy barrier [22]. One expects a crossover 
from thermal activation to quantum tunneling to occur when 
[23], in (55) in place of thermal activation we use the 

quantum tunneling: k B T =h/T GL . A vortex moves 
from x = — > W , during the time T GL . 

We estimate the total interaction energy interaction with 
the neighbourhood vortexes or with the one giant-vortex, 
Fig. lb, of others nucleons from the nuclide nucleus, during 
the time T GL along the vortex path by matching (56), (57) 
and (59) as : 

G« = = 

<I> 2w , sx&„h 

^—k(\-k 2 ) °; = (63) 
c 8x22, 

= a -?—^H — 

2 c I 



where the ratio 



■ = z was chosen as a variable in (55), 



t 2 

through (60). In fact, this is the work done by the Lorentz 
force on the vortex path of the length w . 

Now, we proceed to application to some radionuclides 
which decay beta, and begin with the neutron. 

The lifetime of the free neutron is a basic physical quantity, 
which is relevant in a variety of different fields of particle and 
astrophysics. Being directly related to the weak interaction 
characteristics, it plays a vital role in the determination of the 
basic parameters like coupling constants or quark mixing 



angles as well as for all cross sections related to weak p — n 
interaction. From the most precise measurement within this 
class of experiments, we got the result of r n = 886. 3[s] ■ 

From "Chart of Nuclides 2010", we can get the result of 



Mo ( P~ , T l/2 = 65h, Q = 1 .356M<? V ); 



T l/2 =4Sh 



Q p = 0.687 MeV and 
y - ray = 514fey(0.46%) -> T l/2 = 10.756^ 
r - ray = \5\KeV(95%) -> T 1/2 = 4.48/z 
'l 7 5 Cs -> "iBa + e-+v , Q p =\ .\15MeV 
y-ray = 661foV(85%) -> T l/2 = 30.Syr ,etc 

Now numerically, with these data, we have the result of 
t gl =1.5e-24[j] , 7 = 2560 , v h =\51 , and with 

W = rX as variable, the evolutions of dark count rate R , 

and of R n , for different isotopes are given in Fig. 3. 

Here, the fraction of the bias current to critical current 
l/l Q as used in Y(z) from (60) was deduced separately for 

'55 Cs as that of the plateau zone in Fig. 4, respectively 

of // jU 2 I = 0.005, by using the condition R v It 1 ^ 2 = 1 . 

We can observe that this value corresponds with the expected 
value from quarks transformation 

of =\(e ± )* v vonex * \.6e-\9r [Am] ^ 124[ A] , where the 

vortex crossing velocity is V vonx = X/t cl , and r Q -K shell 
radius. Therefore, the bias current, which is perpendicularly 
on the monopoles current, is / = 0.005 /// = 150[A] , 

where, /j 2 — 0. 1 1 1 , and the monopoles current 
is I = 3.e5[A] as given by (57), and I c =1.6e4[A] from 
(58). From (62), the results will be Q hM = 2.92e-13[7] or 
near equally with of the almost of nuclides, for example, 
for 10() Mo, Q fl = 1.356MeV H> 2.17e - 13[7] . 

Now, in case of quantum tunnelling, the transmission 

Vol 

coefficient [20] is T = e " =l.e-66, which is too small, 
meaning that a vortex crossing from one strip edge to the 
opposite one induces a phase slip, which is the only viable 
mechanis m for dark counts (beta decay). 

Thus, it has established a logarithmic equation of the 
J3 decay rate which resulting a straight line as a function of 

the barrier width (w = rX) for every nuclide, as shown in 

Fig. 5. It decreases in case of long lived nuclides, like 60 Fe . 
Therefore, this evolution is a decisive validation test of 
entirely model. The factor r describes the interactions of the 
nucleons inside of the nucleus, being a complex function 

of Z,N,A,Q bmd . 
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w-bamer width=rlambda 
-Rcbk Qjnd -"-RChn 



Fig. 4 The evolution of dark count ( f} decay) rate as function of barrier 

width 




Fig. 5 The evolution of the dark count rate ( /? ) decay as a function of the 
bias current 

B. Who is in Fact the Parameter w = rX? 

In order to answer to this question, in the following we 
will look to the vortexes interactions [24] vis-a-vis of the 
nucleus substructure models [25 — 33] . 

Firstly, from [24], for vortex-vortex interaction (V-V) at 
large separation, using the asymptotic form of the modified 
Bess el functions, it is: 



Q(d -> oo) = yd-' 12 (Se d - V/^"" ) 



(64) 



where d(A) = rX is the separation distance; jiu = 
and y — -12.503; 5 = 0.3127 are fitting parameters. 

Furthermore, in [24], it shows that the interaction between 
an antivortex and a giant vortex is always attractive as well. 

In the absence of lateral confinement, a giant vortex is a 
stable (unstable) state in type-I (type-II) systems and can 
interact as such with other vortices, and this motivates us to 
investigate the interaction force between a vortex and a giant 
vortex. 

Also in [24], it is suggested that the fitting function given 
by the following equation can not only be used forthe V —Vs , 
but also for the V-GV interaction force. 



a r M) = i, — r Ue'-f^) 

l + T) 2 d 2 



(65) 



where jj, (i = 1, 2, 3, 4) are four fitting parameters . 

The values of the four fitting parameters are given in 
Table II of [24], for ju from 0.2-2.5 . In our 

case, ju — yflu: = 1.4 results in 

77, =1.237e-02;?7 2 = 1.55e-03;?7 3 = 4.248; t] 4 =1.183 

When two vortices or a vortex and a giant vortex are 
brought close to each other, they merge forming a single giant 

vortex state with vorticity n = n l + n 2 , and in the limit of 

small separation, V —V or V — GV forces will be very weak. 
Conversely, a vortex and an antivortex attract and annihilate 
both in type I and type II superconductors. As follows, from 
theory [24], the V-AV interaction is always attractive. 

However, at some critical V — AV 



or 



separation, d E = 0.337 +3 1.249(1 + 10.264//)"° 6855 
d E — 4.96A becomes unstable as the V — AV separation is 
reduced at d < d A , d A = 0.337 +12.222(1 + 2.461//)" 7 " or 

d A =3.95/1 . For molecular dynamics studies of the 
V — AV motion, one should consider the critical separation 
d . as the separation where the V — AV pair annihilates. 

The fit [24] of the numerically obtained V-AV 
interaction force for d > d E , yields a single expression: 



(66) 



Q = -(2.879 + .415/r 3 '^K^cO - 
(0.2258 - 1 .044e 1 866 " )//K, (jud) 
where the approximation of the Bessel function is 

71 



K (d) 



exp(-d) i which is expected to provide an 



accurate description of the V-AV interaction force, at 
separations d > d E , for any value of ju . 

Now, atomic nuclei is known to exhibit changes of their 
energy levels and electromagnetic transition rates among them 
when the number of protons and/ or neutrons is modified, 
resulting in shape phase transitions from one kind of 
collective behaviour to another. These transitions are not 
phase transitions of the usual thermodynamic type; they are 
quantum phase transitions QPT [25] (initially called ground 
state phase transitions). 

A main conclusion from [25] is that the effect of the 
fermionic impurity ( e ± as in case of /3 ± decay) is to wash out 
the phase transition for states with quantum numbers: 
K = 1/2,3/2,5/2 and to enhance it for states 
with k = 7/2, 9/2, ll/2 ■ 

An important property of atomic nuclei is that they 
provide experimental evidence for shapegpjj., in particular, 

of the spherical to axially-deformed transition (U(5) - SU(3)) ■ 
One of three signatures has been used to experimentally verify 
the occurrence of shape phase transitions in nuclei, namely: 
the behavior of the gap between the ground state and the first 

excited + state. 

The nucleus 152 Sm , with n = 90 and z = 62 , lies 
intermediate between nuclei of known spherical shape and 
well-deformed axially-symmetric rotor structure. A sudden 
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change in deformation occurs at N » 88-90 for the Sm and 
neighbouring isotopic chains, see Figure 6. 

New data obtained in [31] constrains the description of 
this nucleus within the IBM to parameter values near the 
critical point of the transition from oscillator to rotor structure. 

The performed IBA calculations in [31] for the entire 
region N « 90 with the simple IBA Hamiltonian which 
involves only one parameter, £ which depends only on the 

neutron number N , for all the isotopic chains, are given in 
Figure 6, where is shown that the evolution of some basic 
observables in Nd — Sm is very well reproduced by these 

calculations. The energy of the intrinsic excitation + has a 
minimum at the phase transition point. The phase/shape 
transition is mirrored in the calculations and, as noted, the 
IBA parameter £ plays the role of a control parameter. 

The total energy surface [31] corresponding to the Sm 
isotopes obtained from the IBA with increasing neutron 
number changes the location of the deformation minimum, 
from (5 — to finite p when the neutron number increases 

from 88 to 92. 

Notice that the bosons (Z(62-50)/2 = 6;N(90-82)/2 = 4) 

pairing appears for l ^Sm 90 (stable) <- EC <- l ^Eu w , 

and ™Sm** (stable) <- yT <- ^Pm 89 , until a free proton 

appears for ^Sm" + B~ -> l ^Eu && (stable) . and 

'115m 9 1 ^>/3-^> l §Eu 92 (stable) ; l gSm 93 -> JB~ -> 'ff £w 92 > 
the same situation for 



Mi 8 



>fi 



Pm 81 



0- -> l gSm u (stabk)>efc- 



Now, in terms of our model defined interactions, it means 
that a proton is composed of 2 vortexes (y) and 1 antivortex 

(AV) of different spin, a neutron is a Giant-vortex (GV) ■ 
Therefore, for TV < 90 , we have an 
interaction (y _ y) _ (y _ AV) > an d above of 
(y -AV)-(V -GV) type, respectively. Based on these 
rationales we have obtained the values of£(0 + ), as shown in 
Fig. 6, where, we show distinctly the values of interactions 
(y -GV);(V -V);(AV -V) as calculated with 
equations (64-66) • 
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Fig. 6 The evolution of empirical values of E(0* 2 ) m me Nd — Sm isotopic 
chains with M > 82 compared with the IBAresults asextracted from [31], 
and of as calculated by the present model 



In both models the phase/shape changes , in our model, 
the phase slip which conducts to a dark count accompanied by 

a W boson energy depends on a single parameter w = r • A , 

till in IBA model are functions of a control parameter £ and a 

quadrupole deformation p as an order parameter, 
respectively, see Fig. 6. 

Our model has the advantage to explain the cause of phase 
change and its associated probability, as expressed by the 
inherently dark counts rate or the decay half-time Ty 2 ■ 

IV. CONCLUSIONS 
A model was developed in which a photonuclear reaction 
is viewed as an incident photon creating superconducting hot 
spot (hot belt) across nucleon from the composition of a 
unstable nuclide (radioisotope), followed by a thermally 
induced vortex crossing, which turns superconducting hot belt 
into the normal state (vacuum) resulting in a vortex assisted 
photon beta decay. Because this model, firstly developed for 
superconducting nanowire single -photon detector (SNSPD), 
requires data on energies of vortex and on currents from 
inside of the nucleon, an analog of a superconductor model 
was developed for the nucleon. 

Thus, a model G-L was revisited in order to calculate the 
Lorenz force, the current, and the energies of the Abrikosov 
vortex lines inside of the nucleon. It was found that these 
energies correspond to subatomic particles, W , Z , Higgs 
bosons, and of meson 7T . Therefore, the nucleon can be seen 
as a triangular lattice with three pairs of quarks -antiquarks in 
the tips of the triangle. These axial vortexes (filaments) 
interacting in the lateral plane. 

All of the keys of this model, as being analogous to a 
superconductor clinging very well, starting with the use of 
Maxwell's equation with monopoles, further mass, charge, 
and the number of monopoles (density), which define the 
penetration depth and the coherence length could be 
considered as new fundamentals constants fornucleons. Also, 
we can say that, because no free quarks were detected, the 
same is true for monopoles, for both of which are confined 
together and in full in nucleons, when the temperature of the 

universe has reached 2* 10" " K . 

Finally the model provides counts (decay) rates ( R ) 
comparable with the experiments, and allows the estimation 
of intensity R j J y photons j ' s} , and of threshold value of 

photons energy ( Me V ) , in order to achieve R pr j R h = 1 , or 

(100%) efficiency of beta decay type radionuclides in the 
frame of the nuclear transmutation in the reducing of 
radioactive waste, together with a decay enhancement 
method. As a validation test, this is the application of this 
model to the mechanism of natural beta decay viewed as dark 
count rate, which was discussed in comparison with the 
results of nuclide substructure model IBA. 
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